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Resumen 

En un cauce natural o artificial es usual la presencia de vegetación en el 
fondo y en la pared. Los efectos que produce la presencia de ésta es una 
alteración en el campo de flujo. La magnitud de las alteraciones en el 
campo del flujo depende de la extensión, altura, densidad y flexibilidad 
de la vegetación. Para analizar a detalle las estructuras turbulentas, en 
este trabajo se realizó un estudio experimental en laboratorio en un canal 
de sección rectangular con un tramo de vegetación sumergida, flexible y 
artificial. La región de medición se ubicó aguas arriba de la zona de 
vegetación hasta el final de esta misma, además se determinó la asimetría 
que puede inducir la presencia de la pared del canal. Las mediciones se 
realizaron con equipo de velocimetría acústica de efecto Doppler (ADV). 
Los resultados obtenidos identifican el desarrollo de la capa de mezcla a 
partir de los patrones de velocidad media y los esfuerzos de Reynolds. En 
el caso de los perfiles de velocidad media sobre la capa de vegetación fue 
posible identificar un punto de inflexión donde se identificaron 
inestabilidades tipo Kelvin-Helmholtz. Con el fin de analizar los procesos 
en la capa de mezcla se propone un modelo fenomenológico basado en el 
análisis de la curvatura de las velocidades promediadas. Además, se 
determinó la aportación en el transporte de momento de las estructuras 


turbulentas con un análisis de cuadrantes con las velocidades fluctuantes. 
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de mezcla, estructuras turbulentas. 


Abstract 

In a natural or artificial channel, the presence of vegetation at the bottom 
and on the wall is usual. The effects produced by the presence of this is 
an alteration of the flow field. The magnitude of the changes in the flow 
field depends on the extension, height, density, and flexibility of the 
vegetation. To analyze the turbulent structures in detail in this work, an 
experimental study was carried out in a laboratory, in a rectangular 
channel with a section of submerged, flexible and artificial vegetation. The 
measurement region was located upstream of the vegetation zone until 
the end of it, in addition, the asymmetry that the presence of the channel 
wall can induce was determined. The measurements were performed with 
an Acoustic Doppler Velocimeter (ADV) instrument. The results allow to 
identify the development of the mixing layer from the average velocity 
patterns and the Reynolds stress. In the case of the average velocity 
profiles over the vegetation layer, it was possible to identify an inflection 
point where Kelvin-Helmholtz instabilities type were identified. In order 
to analyze the processes in the mixing layer, a phenomenological model 
is proposed based on the analysis of the curvature of the averaged 
velocities. In addition, the contribution to the moment transport rate by 
the identified turbulent structures were determined with a quadrant 


analysis of fluctuating velocities. 
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Introducción 


En un cauce artificial, la interacción del flujo con la vegetación produce 
una zona de estelas y se generan regiones de producción de turbulencia 
que alteran los patrones de flujo medio, con una tasa tridimensional muy 
alta. 

En el caso de un cauce natural, la presencia de vegetación 
sumergida es normal, pero los cambios en su morfología están 
relacionados con su presencia e interfieren en función de su altura, 
densidad y flexibilidad. La vegetación genera cambios en las estructuras 
turbulentas en sentido vertical y transversal, además de flujos 
secundarios importantes. Para evaluar la magnitud de estos cambios es 


recomendable estimar patrones de flujos medios y esfuerzos de Reynolds, 
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como una condición de caracterización de baja frecuencia, pero además 
se deben evaluar estadísticamente las velocidades fluctuantes. 

En esta parte del documento se describen las condiciones de 
sumergencia, flexión y densidad de la vegetación. La condición de 
sumergencia, según Palmer (1945), la establece como la relación entre la 
altura de la vegetación respecto al tirante del cauce. La composición 
interna de ésta es muy importante, ya que influye de forma muy 
específica en la modificación del fondo del canal. Asimismo, un parámetro 
importante a considerar debe ser la relación entre la respuesta de la 
flexión de los elementos de vegetación en relación con la cinética del 
campo de flujo. 

La flexión de la vegetación (J) es la respuesta dinámica de la 
estructura de los tallos de la vegetación a la acción de la velocidad 
principal del flujo (Okamoto 8 Nezu, 2010; Ghani, et al. , 2013; Jeon, 
2015), la cual se encuentra en función de la velocidad principal y se puede 
clasificar de acuerdo con los resultados experimentales obtenidos por 
Carollo, et al. (2005), y Okamoto 8: Nezu esta clasificación se enlista a 


continuación: 


I. Rígido. La vegetación se mantiene en su posición vertical y no 
presenta oscilaciones en la copa. 
II. Oscilación suave. La vegetación tiene una deformación suave y 


tiene una oscilación suave no organizada. 
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III. Monami. La vegetación tiene una deformación significativa, y se 
observa un movimiento oscilatorio coherente y armónico. 
IV. Propensa. La vegetación tiene un abatimiento de todos sus 


elementos y las oscilaciones son pequeñas con alta frecuencia. 


Adicional al comportamiento a la flexión, la vegetación se puede 
clasificar respecto a su ubicación y distribución, que puede ser en la pared 
o fondo del canal exclusivamente, o en ambos lugares. También es 
importante establecer la densidad de la vegetación, que puede inducir 
diferentes comportamientos para un mismo tipo de relación entre la 
velocidad y la flexión (Poggi, et a/. 2004; Nezu 8 Sanjou, 2008). 

La relación entre la altura de la vegetación y el nivel de agua, 


definida por Nepf 8 Vivoni (2000), se clasifica en tres tipos: 


I. Flujos bajos: el flujo que pasa a través de los tallos y el follaje son 
deformados. 
IT. Flujos intermedios: se tiene una sumergencia completa de la 
vegetación, pero afecta al comportamiento del campo de flujo. 
III. — Flujos altos: el campo de flujo no se ve afectado por la presencia 
de la vegetación; para esta condición, la vegetación se considera 


como una rugosidad de fondo uniforme. 


El estudio de flujos intermedios se ha convertido en un tema 
importante a analizar debido a los efectos que se generan en el flujo 
medio y turbulento, ya que modifican los perfiles de velocidad media, 


generando estructuras turbulentas coherentes cerca de la copa de la 
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vegetación; estos vórtices son principalmente de gran escala, provocando 
una zona de intercambio de masa y momento entre la copa de vegetación 
y zona exterior, es decir, sobre la vegetación, una condición semejante a 
una capa límite. En la práctica de la ingeniería, la vegetación sumergida 
se analiza como una condición de resistencia uniforme, similar a un flujo 
alto. Entonces, para evaluar su efecto, se aplica el concepto usual de la 
resistencia del flujo por efecto de la vegetación sumergida y sus efectos 
en los parámetros de rugosidad, obteniendo valores generales de 
coeficiente de rugosidad de Manning en función de la relación de 
sumergencia (Járvelá, 2004; Wilson, 2007). 

En este trabajo se realizará un estudio más localizado de los efectos 
de la vegetación sumergida, considerando que el desarrollo del campo de 
flujo sobre la copa de la vegetación es similar a la condición de dos 
campos de flujo a diferente velocidad, con una capa de interacción entre 
éstas y que produce pequeñas perturbaciones, que inducen una 
inestabilidad del flujo del tipo Kelvin-Helmoholtz (KH). Por lo anterior se 
generan vórtices de gran escala, como muestra Finnigan (2000). Para 
analizar con mayor detalle la generación de las estructuras turbulentas, 
en principio se puede decir que se tiene una capa límite y esta se 
desarrolla longitudinalmente, pero Ghisalberti € Nepf (2002) demuestran 
que los flujos sobre vegetación sumergida pueden analizarse como una 
capa de mezcla, más que como una capa límite. El perfil de velocidad de 
una capa de mezcla se aproxima a la forma de una curva tipo tangente 


hiperbólica (Ho 8 Huerre, 1984), como se muestra en la Figura 1. 
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Figura 1. Perfil de velocidad en una capa de mezcla (Ghisalberti 8 Nepf, 
2002). 


De acuerdo con Poggi et al. (2004), en la zona de vegetación se 
pueden identificar en tres regiones, según el tamaño de los vórtices, como 
se muestra en la Figura 2. La primera región está dominada 
principalmente por vórtices de pequeña escala (vórtices de von Kármán). 
La segunda región, que está cerca de la copa de la vegetación, está 


dominada por los vórtices de Kelvin-Helmholtz generados por la 
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inestabilidad en el punto de inflexión. La tercera región es similar a una 


zona de desarrollo tipo capa límite. 


Fo REGIÓN IN 
' $ 5 REGIÓN Il 
OL 5 


REGIÓN I 


7 


Figura 2. Tres escalas de vórtices en un canal con vegetación 


sumergida (Poggi et al., 2004). 


Sukhodolov € Sukhdolova (2006) realizaron mediciones en campo 
para verificar la aplicabilidad de la analogía de la capa de mezcla en un 
cauce natural sobre un parche de Sagittaria sagittifolia implantados con 
un cierto arreglo variando su densidad. Los resultados confirman la 
eficacia de la analogía de la capa de mezcla. Sukhodolov € Sukhdolova 
(2006), y Maltese et al. (2007) se enfocaron en el patrón espacial de la 


estructura turbulenta que se desarrolla sobre la vegetación sumergida. 
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Además, aplicaron el análisis de cuadrantes, revelando que la estructura 
dominante es la expulsión sobre la copa de vegetación. 

Ghisalberti 8 Nepf (2002) revelaron que los flujos con vegetación 
sumergida se pueden modelar como una capa de mezcla en lugar de una 
zona de capa límite debido a que los perfiles de velocidad que se 
presentan en un canal con vegetación sumergida se aproximan a un perfil 
de tipo tangente hiperbólico (Figura 1). Además, a partir de los perfiles 
de velocidad media identifican un punto de inflexión, lo cual implica que 
el flujo es susceptible a las inestabilidades de Kelvin-Helmholtz. Nikora et 
al. (2008) examinaron el efecto de la vegetación sumergida en la 
resistencia hidráulica en una variedad de patrones de parches de 
vegetación. Okamoto, et al. (2012) realizaron experimentos con inyección 
continua de tinta para evaluar la estructura del transporte de masa en un 
flujo de canal abierto con modelos rígidos de vegetación y una 
configuración de tiras rígidas variando la densidad, y obtuvieron los 
efectos de los vórtices coherentes en la difusión turbulenta vertical. 
Okamoto 8 Nezu (2013) examinaron la transición de una capa límite 
desarrollada aguas arriba del parche de vegetación a un flujo de tipo de 
capa de mezcla desarrollada sobre la vegetación; además desarrollaron 
un modelo fenomenológico para el proceso de desarrollo en el flujo con 
vegetación sumergida. En estos experimentos utilizaron la técnica PIV 
(Particle Image Velocimetry) y los elementos de vegetación estaban 


compuestos por tiras rígidas. 
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Como ya se mencionó, las características de flujo medio y 
turbulencia en flujos de canal abierto con vegetación sumergida han 
recibido mucha atención. Sin embargo, aún no se dispone de información 
detallada de la evolución de la capa de mezcla y las estructuras 
turbulentas coherentes generadas por la presencia de la vegetación. Por 
esto, en el presente estudio se considera el efecto de un parche de 
vegetación flexible de longitud finita y un flujo intermedio para identificar 
las estructuras turbulentas, al igual que su evolución longitudinal y 
transversal a partir de un modelo fenomenológico. La zona de análisis 
inicia aguas arriba del parche hasta el límite del parche, pero se realizaron 
mediciones en secciones transversales y verticalmente se consideró la 


zona por debajo de la copa de vegetación hasta el flujo exterior. 


Materiales y métodos 


La estación experimental es un canal con pendiente fija S, de sección 
rectangular, con una longitud de 10 m, ancho de 0.60 m y con una altura 
de 0.60 m. Las paredes del canal son de vidrio y fondo de acero. La 


representación de los elementos de vegetación es de un material plástico 
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flexible. Las condiciones hidráulicas del estudio experimental se muestran 


en la Tabla 1. 


Tabla 1. Condiciones hidráulicas donde Q es el gasto; U = Q/A es la 

velocidad media del flujo; Sy es la pendiente de la plantilla del canal; 

Re = UR;,/v es el número de Reynolds; Fr = U/y/gR, es el número de 
Froude; R, = A/P es el radio hidráulico; A es el área de la sección 
transversal; P es el perímetro de la sección transversal; v es la 


viscosidad cinemática, y g es la aceleración de la gravedad. 


h(cm) | h,(cm) | h/h, | Q(s%) | U(cms”) S, Re Fr 


33 11 3 41.20 20 0.005 | 32698 | 0.11 


El sistema de referencia de la toma de muestras de la estación 
experimental se muestra en la Figura 3, donde el sistema coordenado es 
la tríada de valores de x;,, para ¡=1,2,3, y se denominará coordenada 
longitudinal, vertical y transversal, respectivamente. El dominio del 


espacio de medición en metros es M(x;) € [-0.5, 3.7] x [0.0, 0.33] x [0.0, 0.3]. 
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Figura 3. Estación experimental. 


Para la toma de las velocidades instantáneas en las tres 
componentes (u;) se utilizó un dispositivo de velocimetría acústica de 
efecto Doppler (ADV, Acoustic Doppler Velocimeter) a una frecuencia de 
muestreo de 100 Hz durante un periodo T = 60 s. El dispositivo ADV es 
de la marca Nortek) o también llamado Vectrino Profiler”" (Vectrino IT). 
Para el proceso de los datos experimentales se consideró el modelo de 
descomposición de Reynolds u; = U; +u;, considerando los promediados 
temporales de las velocidades en cada componente como U;, y las 
velocidades fluctuantes como u;. 

La configuración de la vegetación es un parche con una longitud de 
3.70 m y el arreglo de elementos flexibles es lineal, con una separación 


entre elementos en sentido longitudinal de 1, =0.025m y en sentido 
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transversal de !l,, =0.025m, tal como se muestra en la Figura 4 y Figura 
5, respectivamente. El inicio del parche de vegetación se considera en 


x1, =0. 


X2 
Flujo 


Zona de fondo liso a 


h =0.33 m 


—=—).50 SA A o X1 


Figura 4. Zona de vegetación. 
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X3 


X1 


Figura 5. Patrón de los elementos de vegetación. 


Los puntos de muestreo fueron del centro de canal hacia la margen 
derecha en cinco secciones, indicadas con las letras A, B, C, D, E (Figura 
6); en sentido longitudinal se analizaron 43 puntos y 19 en sentido vertical 
(Figura 7a). Los puntos de medición en la vertical se muestran en la Figura 
7b. 
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Figura 6. Mediciones en sentido transversal. 


a) 


Figura 7. Puntos de medición con el equipo ADV, con un total de 4 085 
puntos: a) puntos de muestreo en el plano (x,,x3) con una separación 
de 0.10 m en sentido x, y 0.06 m en sentido x3; b) puntos de muestreo 
en el plano (x;,, x,), con una separación de 0.015 m en sentido x,. En la 
zona de vegetación se tienen siete puntos en el interior y 12 encima de 


la vegetación. 


42 
2023, Instituto Mexicano de Tecnología Tecnología y ciencias del agua, ISSN 2007-2422, 
del Agua. Open Access bajo la licencia CC BY-NC-SA 4.0 14(2), 27-82. DOI: 10.24850/j-tyca-14-02-02 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


a 0) Check for updates 
OPEN ACCESS 
Tecnología y 


CienciaszAgua 


De acuerdo con Vargas et al. (2015), la vegetación real puede ser 
simulada con barras cilíndricas. Con esto, se propuso utilizar barras 
cilíndricas para simular la vegetación con las siguientes características 
físicas: altura de hp = 0.11 m y diámetro de / = 7.3 mm. El material de los 
elementos de vegetación es silicón. 

Para la caracterización de la vegetación se utilizaron los parámetros 
de densidad de vegetación 2A;, rigidez a la flexión J y la relación de 
sumergencia h/h,. La densidad de vegetación adimensional 2; se 


determinó con la expresión propuesta por Okamoto 8: Nezu (2013): 


En E" Af ¡hp 
y == (1) 


Donde 4 es el área frontal total de los elementos de vegetación y 
A, es el área ocupada por la vegetación. En el caso de la rigidez a la 
flexión (J) de los elementos de vegetación se utilizó la ecuación propuesta 


por Wilson, et al. (2003), que es: 
J=El === (2) 


Donde E es el módulo de elasticidad e / es el segundo momento del 


área (también llamado como momento de inercia). En el lado derecho de 
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la expresión se tiene que F es la relación entre la fuerza; ó,,, es la deflexión 


del elemento; L, la longitud del elemento de vegetación (Figura 8). 


Entonces, se obtiene un valor de J = 3.42x10"? Nm? y un valor A, = 1.37. 


X2 


X1 


Figura 8. Diagrama del elemento de vegetación con deflexión. 


Para este estudio se simuló un tipo de vegetación real llamado 
Lespedeza sericea, con una longitud de h, = 10 — 20 cm y con un J = 5x107* 
a 1x107? Nm? (Kouwen € Li, 1980; Nepf 8 Vivoni, 2000). Los parámetros 
utilizados en este experimento, tales como la relación de sumergencia, la 


densidad, la rigidez a la flexión y el patrón de los elementos de vegetación 
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fue con base en varios experimentos de laboratorio realizados por 
diversos autores (Palmer, 1945; Wilson et al., 2003; Nezu € Sanjou, 
2008; Okamoto € Nezu, 2013). Los patrones lineales y escalonados son 
los más utilizados en laboratorio debido a que tienen una mejor 
representación del efecto de la vegetación en un canal con flujo 


intermedio (Vargas-Luna et al., 2015). 


Resultados 


De acuerdo con lo anterior, la condición de flujo que se está analizando 
es un flujo intermedio. Además, para una mejor comprensión de los 
resultados se consideró un espacio adimensional del dominio de muestro 
O(¿; = x¡/hp) € [-4.54, 33.63] x [0.0, 3.0] x [0.0,2.72]; la copa de la vegetación 
se tiene para ¿, = 1, y el inicio de la vegetación para ¿, = 0. Los valores 


para las secciones en sentido transversal se muestran en la Tabla 2. 
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Tabla 2. Valores de ¿3 para los puntos de medición en el sentido 


transversal. 
Sección A B C D E 
¿3 =x3/h, 2.72 2.18 1.63 1.09 0.54 


Características de flujo medio 


De las mediciones obtenidas, primero se analizaron los campos de flujos 
medios. En la Figura 9 se muestran los contornos de velocidad 
promediada temporal en el sentido principal U,. La velocidad promediada 


en los puntos de muestreo (¿;) se determinó con la siguiente expresión: 


UE) =5S0 uE, man (3) 
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$3 = 2.72 Es = 2.18 


U¡ (m/s) 


Figura 9. Contornos de velocidad media en sentido principal U,, donde 


la línea discontinua representa la copa de la vegetación £, = 1. 


Por la presencia de la vegetación se observa una desaceleración de 
la velocidad principal cerca de la copa de vegetación ¿,=1, que se 
compensa con un incremento de ésta en la zona exterior. Esta distorsión 


del campo de flujo origina un punto de inflexión justo en la copa de 
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vegetación ¿,=1, lo anterior genera una inestabilidad de Kelvin- 
Helmholtz y por consiguiente se presentan vórtices de diferentes escalas; 
esta observación se asemeja a resultados presentados por Nepf 8 Vivoni 
(2000). 

En la Figura 10 se muestran los perfiles de velocidad media para 
todas las secciones transversales, observando un perfil tipo hiperbólico. 
Dichas características del campo de flujo son similares a las que se 
presentan en la vegetación terrestre, presentada por Finnigan (2000). La 
inestabilidad se mantiene en el sentido transversal, pero disminuye el 


espesor de la capa de mezcla por la presencia de la pared. 
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E, = 1181 E, = 17.27 E, = 22.72 E, = 28.18 E, = 33.63 


$1 = 4.54 E, =0.90 E, =6.36 


Sección transversal 
—+— Ez = 1.63 

—— Ey = 2.72 —t— Ez = 1.09 

—A— Ez = 2.18 —*— ¿Ez = 0.54 


Js 
0 0.2 04 0 0.2 


U, (m/s) 


Figura 10. Perfiles de velocidad media en sentido principal, U,, donde la 
línea discontinua indica la copa de la vegetación. 


En el perfil de velocidad media en é¿, = -4.54 (Figura 10) en la zona 
pegada a la superficie libre (¿, > 2.5) se observa una aceleración del flujo 


por el efecto de remanso ocasionado por la vegetación. 
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Esfuerzos de Reynolds 


Analizando los esfuerzos de Reynolds —u;u;, de los experimentos (Figura 
11), se puede observar el desarrollo de la capa de mezcla, que es muy 
representativa para la zona central (¿3 = 2.72), pero ésta va disminuyendo 


a medida que se tiene la influencia de las paredes, ¿z < 2.72. 
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$3 = 2.72 ¿2 = 2.18 


0 5 10 15 20 25 30 E ) 


Figura 11. Contornos de esfuerzos de Reynolds —u;u; , donde la línea 


discontinua representa la copa de la vegetación ¿, = 1. 


Nepf 8: Vivoni (2000), Nezu 8: Sanjou (2008), definen el tamaño de 
la capa de mezcla a partir de los perfiles de esfuerzos de Reynolds —u;u; , 
donde el límite superior e inferior están definidos por el 10 % del esfuerzo 


máximo [uu es decir, como se muestra en la siguiente expresión: 
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_ 1 o(|-uuz),..) (4) 


Donde y, es el límite superior y y, es el límite inferior de la capa de 
mezcla. Los límites superior e inferior de la capa de mezcla para ¿z = 2.72 
y ¿3 = 0.54 se muestran en los perfiles de esfuerzos de Reynolds (Figura 
12). De acuerdo con el perfil de los esfuerzos de Reynolds para la posición 
longitudinal ¿, = 33.63 y en la sección transversal ¿z = 2.72, se tiene un 
espesor de capa de mezcla de t,»¡, = 0.19 m y para ¿ = 0.54 de t,,,, = 0.15 m. 
Con esto se observa que al centro del canal se tiene un tamaño de capa 


de mezcla mayor al que se tiene pegado a la pared. 
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E, = 4.54 E, = 0.90 E, = 6.36 E, = 1181 E, =17.27 E, = 22.72 E, = 28.18 E, = 33.63 


Uma 


Sección transversal 
—— [q = 1.63 

—e— E, = 2.72 —+— Ey = 1.09 

—a— E = 2,18 —+— Ez = 0,54 


Ly Ly 
10 15 10 15 
x 10% «10% 


uu, (m?/s?) 


Figura 12. Perfiles de esfuerzos de Reynolds —uu), donde las líneas 
discontinuas color azul y magenta representan la variación espacial en 
sentido longitudinal del espesor de la capa de mezcla en la sección 
central ¿¿ = 2.72 y pegada a la pared ¿3 = 0.54, respectivamente, y la 


línea discontinua color negro representa la copa de la vegetación ¿, = 1. 
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Capa de mezcla 


Según Raupach, et al. (1996), y Ghisalberti 8 Nepf (2002), el perfil de 


velocidad de una capa de mezcla se puede analizar con la siguiente 


expresión: 
Ui-Um _ 1 (E2) 
a. > ¿tanh EA (5) 


Para una mejor comprensión en la Figura 1 se muestran las 
variables de la capa de mezcla. Donde U,, = U,(y,) + U,(y,)/2; AU, es la 
diferencia de velocidad para cada sección transversal: AU, =U,(y,) — 
0.0); U0,) y 0,07, ) indican la velocidad promediada en el límite superior 
e inferior de la capa de mezcla, respectivamente; y, y y, es la distancia 
del fondo del canal al valor del esfuerzo de Reynolds con valor cero por 
debajo y encima del valor máximo de esfuerzo de Reynolds, 
respectivamente, y = (y, + y,)/2. Este análisis se desarrolla para cada 
perfil de velocidad en las secciones longitudinales y transversales. El 


espesor de momentum 8 se define con la siguiente expresión: 


o rl] 6 
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En la Figura 13 se hace una comparación entre el modelo teórico 
para evaluar los perfiles de velocidad de tangente hiperbólica de una capa 


de mezcla pura y los datos experimentales. 
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$1 = 0.90 E, = 6.36 


pa mi z 4 L pa n z 
1 03 06 04 02 0 02 04 06 08 1 


$1 = 11.81 E, = 17.27 


xy 


A 40.5 0 0.5 1 


J 5 a 
1 A 0.5 o 0.5 1 


— Tangente hiperbólico 
e ¿=272 AU, 
mn ¿¿=218 
* ¿ =1.63 
% ¿¿=109 
4 £,=0.54 


Figura 13. Comparación de los perfiles de velocidad media con el perfil 
de tangente hiperbólica de una capa de mezcla para cada sección 


transversal para diferentes posiciones en sentido longitudinal. 
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De acuerdo con los resultados, se observa que los perfiles de 
velocidad obtenidos en la experimentación se asemejan a un perfil de 
capa de mezcla, por lo tanto se concluye que se forma una capa de mezcla 


en un flujo intermedio. 


Modelo fenomenológico 


Para evaluar la evolución de la capa de mezcla en sentido longitudinal 
para cada sección de análisis se propone aplicar una división del campo 
de flujo en cuatro zonas, como se describe en la Tabla 3 y se muestra en 


la Figura 14. 
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Tabla 3. División de las zonas de modelo fenomenológico, donde X, es 
el inicio de análisis (zona aguas arriba de la vegetación); X»,, el inicio de 
la zona divergente; Xp, el inicio de la zona de desarrollo; Xzp, el inicio 


de la zona completamente desarrollada, y X, es el fin del campo de 


análisis. 
Zona Condición de flujo Espacio de 

análisis 

I Inicial [Xi Xp] 

II Divergente [Xp Xpp] 

III Desarrollo [Xp Xzp] 

IV Completamente desarrollado [Xzp, X1,] 
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Figura 14. Modelo fenomenológico del desarrollo del flujo a través de 
una zona de vegetación sumergida (Okamoto 8: Nezu, 2013; Huai 8 
Qian, 2015). 


Aguas arriba de la zona de vegetación (¿, < X»,), los perfiles de 
velocidad se asemejan a un perfil encontrado en un flujo canal abierto, 
ya que la vegetación no afecta el flujo. En el inicio de la vegetación (€, = 0) 
existe una alta resistencia al flujo, e inicia un incremento de la velocidad 
por encima de la vegetación y una desaceleración dentro de ésta. Es decir, 


que parte del flujo se desvía al inicio de la zona de vegetación. Zong € 
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Nepf (2010) denominan este fenómeno como flujo divergente. El fin del 
flujo divergente es denotado como Xp. 

Aguas abajo de la zona de flujo divergente (¿, > Xp), los perfiles de 
velocidad tienen un alto gradiente de velocidad, como se muestra en la 
Figura 1. Este gradiente crea un punto de inflexión en ¿, = 1, induciendo 
una inestabilidad tipo Kelvin-Helmholtz, generando vórtices aguas abajo. 
El crecimiento de estos vórtices se observa en los perfiles de distribución 
de los esfuerzos de Reynolds. Esta zona se define como zona de desarrollo 
y el fin de la misma se define como Xp. 

Aguas abajo de la zona de desarrollo (¿, > Xzp) el tamaño de los 
vórtices se mantiene constante y se denomina zona completamente 
desarrollada. 

En la Figura 15 se muestra el efecto de la vegetación en la velocidad 
principal promediada en la vertical U,, en la zona sin y con vegetación. 


Para determinar la velocidad promediada se utiliza la siguiente expresión: 


Ub == 45" U¡ (End, E = ctedag (7) 
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Figura 15. Variación de la velocidad media en sentido principal, aguas 


arriba y sobre la vegetación. 


Observando en la zona aguas arriba de la vegetación, -5<é¿, <0, 
se tiene una velocidad con poca variabilidad, 0U,,,/dx = 0. Para la zona 
central del canal ¿¿ = 1.63, 2.18,2.72, la velocidad del flujo incrementa hasta 
un valor máximo en un rango 0<é, <5, aguas debajo de este valor 


máximo la velocidad disminuye debido a que la capa de mezcla alcanza 
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la superficie libre, como se muestra en la Figura 11, siendo una zona de 
alta turbulencia. Sin embargo, para la zona pegada a la pared ¿;z = 
0.54, 1.09, el incremento de la velocidad se mantiene constante, debido a 
que la capa de mezcla no alcanza la superficie libre, por lo tanto, la zona 
de turbulencia es menor y no se genera la desaceleración del flujo (Figura 
11). 

El punto de inicio del flujo divergente, X»,, y la longitud de la zona 
divergente (Zona IT), Lsp = Xsp — Xp, se determina a partir de la velocidad 
vertical promediada en la vertical, U,, . Aplicando la Ecuación (7) se evalúa 
el flujo divergente, como se muestra en la Figura 16. Se observa un 
incremento rápido de la velocidad vertical en -2.73<€, <0, donde el 
máximo valor se encuentra al inicio de la zona de vegetación (¿, = 0) para 
cada sección. El inicio del flujo divergente, Xy,, se determina cuando la 
velocidad vertical incrementa, es decir, cuando dU,, /dx, > 0, por lo tanto, 
para las cinco secciones Xp, = —2.73. Sin embargo, el inicio de la zona de 
desarrollo, Xp, para cada sección varía a medida que se acerca a la pared. 
Xpp se obtiene hasta el punto donde la velocidad vertical se aproxima a 


cero,U,,, = 0, como se muestra en la Figura 16. En la Tabla 4 se muestra 


la longitud de la Zona Il, Lyp. 
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Figura 16. Variación de la velocidad media en sentido vertical para 


cada sección. 


Tabla 4. Longitud de flujo divergente, L;p. 


Zonaés 2.72 2.18 1.63 1.09 0.54 
XFL 2.73 2.73 2.73 2.73 2.73 
Xp 13.63 12.72 10 7.27 2.72 
Lrp 16.36 15.45 12.73 10 5,45 
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Para visualizar de mejor manera la zona de flujo divergente, en la 


Figura 17 se muestra la zona de mayor velocidad vertical U,. 


$3 = 2.72 Ez = 2.18 


U2¿(m/s) 
0 5 10 


pS 
an 


20 25 30 


ó1 


Figura 17. Contornos de velocidad media en sentido vertical U,, donde 


la línea discontinua representa la copa de la vegetación ¿, = 1. 
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Aguas abajo de la zona de flujo divergente, las velocidades en 
sentido vertical se vuelven pequeñas, observando valores de velocidad 
negativos, U, < 0. Entonces, para ¿, > Xpp, la capa de mezcla crece de 
forma gradual hasta que el tamaño de la capa de mezcla, t,, ,, se mantiene 
constante, es decir, cuando dtm,/dx, = 0, se ubica el inicio de la zona 
completamente desarrollada, Xy». En la Figura 12 se observa la variación 
espacial de la capa de mezcla. Entonces, para cada sección, Xzp se obtiene 
con la variación espacial del tamaño de la capa de mezcla, como se 


muestra en la Figura 18. 
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Figura 18. Variación espacial del tamaño de la capa de mezcla, t.,, 


para cada sección. 


Para comparación, en la Figura 19 se muestra la variación del 
espesor de momentum, 6, que es una medida de crecimiento de la capa 
de corte. Se observa que t,,, y O crecen rápidamente aguas abajo del inicio 


de la zona de vegetación. 
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Figura 19. Variación espacial del espesor de momentum, 6, para cada 


sección. 


Entonces la longitud de la zona de desarrollo se obtiene con Lp = 
Xzp — Xpp. En la Tabla 5 se muestran las longitudes de la Zona III para 


cada sección. 
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Tabla 5. Longitud de zona de desarrollo, Lzp. 


OPEN Que: 


0) Check for updates 


Zonaéz 2.72 2.18 1.63 1.09 0.54 
XD 13.63 12.72 10 7.27 2.72 
Xzp 19.09 20.91 21.82 25.45 27.27 
Lzp 5.45 8.19 11.82 18.18 24.55 


Por último, aguas debajo de la Zona III, €, > Xzp, se encuentra la 


zona completamente desarrollada (Zona IV), L¿p, hasta donde termina la 


zona de vegetación. En la Figura 20 se muestra el modelo fenomenológico 


del desarrollo del flujo para esta prueba, donde se observan 


variaciones de las cuatro zonas en sentido transversal, xz. 
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Figura 20. Modelo fenomenológico, donde la línea discontinua color 


negro indica el inicio de la zona de vegetación, ¿, = 0. 


Se observa que la longitud de la Zona III, L,p, para el centro del 
canal, ¿¿ = 2.72, es menor con respecto a la zona pegada a la pared, ¿z = 
0.54. Esto indica que en el centro del canal existe mayor resistencia, ya 
que en el centro se tiene una desaceleración del flujo alta en comparación 


con la zona pegada a la pared (ver Figura 15). 
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Análisis de cuadrantes 


Se realizó un análisis de cuadrantes para los esfuerzos de Reynolds 
instantáneos —u; (t)us(t) para el análisis de las estructuras turbulentas que 
contribuyen en el transporte de momentum. El principio básico del análisis 
de cuadrantes es descomponer los esfuerzos de Reynolds en cuatro tipos 
de flujos, denominados Q,, Q,, Q3 y Q4, en el plano u; — us, en función de 
los signos de la fluctuación de la velocidad principal u; y vertical us, como 


se muestra en la Figura 21. 


Figura 21. Diagrama de cuadrantes u; — uz. 
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El cuadrante de esfuerzos de Reynolds RS, se define de la siguiente 
manera (Okamoto et al., 2012; Okamoto 8: Nezu, 2013): 


: T 1 
RS; = lim 7/, (ujuz)l; de (8) 


Donde T es la duración de la medición. Si (u;,u,) existen en 
cuadrante Q;, entonces /; = 1; en caso contrario, /; = 0. 
Cada cuadrante de (u;,u,) corresponde a los eventos coherentes 


siguientes: 


Q; (u; > 0,u) > 0): interacción externa, que es el movimiento ascendente 
del flujo de alta velocidad. 

Q) (u, <0,u) > 0): expulsión, que es el movimiento ascendente del flujo 
de baja velocidad. 

Q3 (u, > 0,u) < 0): interacción interna, que es el movimiento descendente 
del flujo de baja velocidad. 

Q, (u, < 0,us < 0): barrido, que es el movimiento descendente del flujo de 


alta velocidad. 


En la Figura 22 se muestra la distribución vertical de los cuadrantes 
de esfuerzos de Reynolds, RS;¡, para las cuatro zonas del modelo 


fenomenológico en ¿¿=2.72. Además, se presenta el plano de las 
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velocidades fluctuantes (u; —u5) para ¿,=1.2, donde se indica el 
porcentaje de contribución de cada cuadrante. En la zona de flujo inicial, 
¿, =-2.72, se observa que los valores son pequeños y el porcentaje de 
aportación de las estructuras turbulentas al transporte de momentum es 


similar, es decir, no se observan estructuras organizadas. 
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Figura 22. Cuadrantes de esfuerzos de Reynolds para la sección central 


del canal ¿3 = 2.72: a) zona inicial (¿, = —-2.72); b) zona divergente 
(€, = 5.45); c) zona desarrollada (¿, = 17.27); d) zona completamente 


desarrollada (¿, = 30). 
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En la zona de flujo divergente, ¿, = 5.45, sobre la zona de vegetación 
(£, > 1) se observa que la estructura dominante es el barrido, RS, (Figura 
22b). 

En la zona de desarrollo, ¿, = 17.27, las estructuras dominantes son 
la expulsión y el barrido por efecto de los vórtices en la capa de mezcla, 
que son los movimientos ascendentes de baja velocidad y descendente 
de alta velocidad, respectivamente Huai et al. (2019). Esto se observa en 
¿, = 1.5, donde se encuentra el valor máximo de esfuerzos de Reynolds 
(Figura 22c). 

En la zona de flujo completamente desarrollada ¿, =30, la 
estructura de expulsión y barrido aumentan un 11 %con respecto a la 
zona de desarrollo. Estos resultados indican que la contribución del 
transporte de momentum entre la zona de vegetación y la zona sobre 
ésta es principalmente por las condiciones de expulsión y barrido. 

Para visualizar las estructuras turbulentas en toda la región de flujo, 
en la Figura 23 se muestran los contornos de los cuadrantes de esfuerzos 
de Reynolds para todas las secciones. Se observa que para todas las 
secciones, las estructuras que mayor contribuyen al transporte de 
momentum son la de estructura de expulsión y barrido. Además, se 


observa que los valores máximos están en la zona de capa de mezcla. 
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Figura 23. Contornos de los cuadrantes de esfuerzos de Reynolds, RS,, 


para cada sección transversal ¿3, evaluado sobre el espacio x, — xz. 


Una vez encontradas las estructuras turbulentas con mayor contribución, 
en la Figura 24 se comparan las contribuciones del movimiento de barrido 


y expulsión utilizando la Ecuación (9): 
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1 ¿S=h RS 
Rsje =58, es, E (9) 


Figura 24. Relación entre las estructuras de barrido y expulsión. 


Aplicando la Ecuación (9) para cada sección se obtienen las 
relaciones entre barrido y expulsión Rs,z, observando una relación Rs,g > 
1 para todas las secciones transversales, por lo tanto, la estructura con 
mayor contribución al transporte de momento vertical se debe al 


movimiento de barrido. 
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Conclusiones 


En este estudio, las mediciones de turbulencia se realizaron en un canal 
con un parche de vegetación de longitud finita para examinar la creación 
y evolución aguas abajo de las estructuras turbulentas en la capa de 
mezcla sobre la vegetación. 

En los perfiles de velocidad media se confirma la presencia de un 
punto de inflexión ocasionado por la inestabilidad de tipo Kelvin-Helmholtz 
similar al análisis teórico cuando se presentan dos flujos potenciales de 
diferente velocidad. Mediante la distribución vertical de los esfuerzos de 
Reynolds se identificó la capa de mezcla, así como la variación del tamaño 
a lo largo del parche de vegetación y en sentido transversal. Se observa 
que el tamaño de la capa de mezcla deja de crecer en la zona 
completamente desarrollada (x, > Xzp) debido a que la capa de mezcla se 
extiende a la superficie libre. Además, el tamaño de la capa de mezcla 
para ¿z = 2.72 (sección A) es de un 2 Y% más grande en comparación con 
¿¿ =0.54 (sección E); esto se observa claramente en los contornos de 


esfuerzos de Reynolds (Figura 11) debido a que el campo de flujo de alta 
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velocidad es más grande al de baja velocidad, como se muestra en la 
Figura 9. 

Con el modelo fenomenológico del flujo se identificaron las cuatro 
zonas del flujo para las cinco secciones de análisis; se obtuvieron las 
longitudes de estas zonas para las secciones, como se muestra en la 
Figura 20. Se observa que para ¿z=0.54 (sección E) requiere mayor 
distancia, para que la capa de mezcla se desarrolle completamente, en 
comparación con ¿; = 2.72 (sección A). Esto se visualiza en la magnitud de 
los esfuerzos de Reynolds. 

Se aplicó un análisis de cuadrantes para identificar las estructuras 
turbulentas que contribuyen al transporte de momentum; para la zona de 
flujo inicial no se observa una estructura organizada; en las zonas de 
desarrollo y completamente desarrollada aparecen estructuras de 
expulsión y barrido cerca de la copa de la vegetación; tales resultados 
indican que la contribución del transporte de momentum entre la zona de 
vegetación y la zona sobre ésta es principalmente por las estructuras de 
expulsión y barrido. 

Con la metodología empleada en este artículo es posible identificar 
la presencia de una capa de mezcla y las características de las estructuras 
turbulentas presentes en ella. Además, con el uso de un modelo 
fenomenológico aplicado a varias secciones se generó un mapa que 
muestra la variación de las cuatro zonas a lo largo de un canal con 
vegetación, lo cual puede emplearse para canales con vegetación con 


características similares a las utilizadas en este experimento, tales como 
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la relación de sumergencia, rigidez a la flexión, densidad adimensional y 


velocidad media principal. 
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